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Nanomechanical detection of the spin Hall effect 
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The spin Hall effect creates a spin current in response to a charge current in a material that has 
strong spin-orbit coupling. The size of the spin Hall effect in many materials is disputed, requiring 
independent measurements of the effect. We develop a novel mechanical method to measure the size 
of the spin Hall effect, relying on the equivalence between spin and angular momentum. The spin 
current carries angular momentum, so the flow of angular momentum will result in a mechanical 
torque on the material. We determine the size and geometry of this torque and demonstrate that 
it can be measured using a nanomechanical device. Our results show that measurement of the 
spin Hall effect in this manner is possible and also opens possibilities for actuating nanomechanical 
systems with spin currents. 


The spin Hall effect [U [2] , which is the generation of 
a spin current in a material due to an applied charge 
current in the presence of strong spin-orbit coupling, has 
been proposed as a novel method of spin manipulation for 
spintronics applications. Spin currents generated by the 
spin Hall effect have been used to excite high- and low- 
frequency magnetic dynamics in nanostructures [SHII], 
and may become useful for future low-power spintronic 
logic and storage devices [12]. The spin Hall effect has 
been observed in a variety of materials with strong spin 
orbit coupling, including semiconductors such as Si and 
GaAs [T3UT6] : graphene with adsorbed impurities El; 
heavy metals with strong spin-orbit coupling such as Pt, 
Ta and W [H [TTJ 1181420] ; and metals doped with large 
spin-orbit-coupled impurities [2T1423] . However, quantifi¬ 
cation of the SHE through fundamental parameters re¬ 
mains a challenge. 

The figure-of-merit for SHE materials is the spin Hall 
angle, &sh^ which is often stated as the proportional¬ 
ity between the magnitude of generated spin current and 
the magnitude of input charge current, \Js\ = | 

El El, where Js is the component of the spin cur¬ 
rent perpendicular to the charge current and Jc is the 
charge current. Measurements and characterization of 
Qsh have, as yet, been limited to methods based on opti¬ 
cal, electrical, and magnetic effects, which require knowl¬ 
edge of Kerr rotation coupling, metallic interfaces, mag¬ 
netic properties, and spin diffusion parameters to quan¬ 
tify &SH accurately [mus]. As such, reported values 
of &SH span orders of magnitude for materials such as 
Pt and Pd jH [271435] . Open questions such as the de¬ 
pendence of Qsh on growth conditions, film thickness, 
impurity level, frequency, and other systematic param¬ 
eters must be addressed both experimentally and theo¬ 
retically. Since the inverse SHE is used to measure spin 
transport due to the spin Seebeck effect and spin pump¬ 
ing, accurate knowledge of the spin Hall angle is impor¬ 
tant for metrology. One intriguing new result implies 
that the spin Hall angle is complex-valued, resulting in 
a phase shift between an applied AC charge current and 
the resulting AC spin current [33] . Additionally, the spin 
diffusion length in spin Hall materials may be correlated 


with the spin Hall angle, indicative of the same physics 
affecting both quantities. How these material parameters 
are related is now an area of intense study. 

Because of the large discrepancies in measured spin 
Hall angle values using electrical techniques, independent 
methods of measuring the spin Hall effect are needed to 
accurately determine the spin Hall angle, the underlying 
mechanisms governing the effect, and the value of spin 
transport parameters such as the spin diffusion length. 
Though very promising for future study, as of yet, no 
quantitative results from metals have emerged from such 
a technique. 

This paper describes a novel method of determining 
SHE properties by measuring the mechanical effects of 
a charge current applied to a SHE material using a 
nanomechanical resonator. This technique relies on the 
equivalence between spin angular momentum and me¬ 
chanical angular momentum, and the fact that one is 
ultimately converted into the other. The equivalence be¬ 
tween spin and mechanical angular momentum has been 
known since the initial experiments of Richardson [36] 
and Einstein and de Haas m, and measurements of 
spin properties through mechanical effects have been per¬ 
formed to study spin injection from a ferromagnet into 
a nonferromagnetic metal [38] [39]. Recent discoveries 
of new spin-based effects, such as the spin Seebeck and 
spin Hall effects, in a variety of materials, from metals 
to semiconductors to insulators, have motivated us to 
study the associated spin physics micromechanically. We 
show that accurate determination of the spin Hall angle 
is possible with this method, which requires no a pri¬ 
ori knowledge of interfacial transparency, spin scattering 
mechanisms, and other unknown quantities. The results 
presented here also offer the possibility of novel excitation 
methods for nanomechanical devices using spin polarized 
currents generated by the SHE and others, including the 
spin Seebeck effect and topological insulators. 

The spin Hall effect is the generation of a pure spin cur¬ 
rent in response to a charge current in a material with 
substantial spin-orbit coupling. The spin current flow 
direction is perpendicular to the charge current and is 
polarized perpendicular to the directions of charge cur- 
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Wire Width: 20 |jm Wire Material: Platinum 
Wire Length: 90 [jm Spin Hall Angle: 0.2 
Wire Thickness: 20 pm 


FIG. 1. (a) Schematic representation of the spin Hall effect. 
Oppositely polarized spins scatter in opposite directions, re¬ 
sulting in a spin accumulation at the wire boundary, (b) 
Graphical representation of the spin accumulation in a cross 
section of the wire. The blue lines and vector plot illustrate 
the accumulation for A^/ comparable to the wire size and the 
red lines represent a much smaller A^/. (c) Deflection (not 
to scale) caused in a platinum wire due to the spin-transport- 
induced torque density. The bottom of the wire is fixed in 
place while the rest may deform freely. 


rent flow and spin current flow. The effect stems from 
spin-orbit and impurity scattering which causes oppo¬ 
sitely polarized spins to scatter in opposite directions, as 
illustrated in Figure [^. This is described using the spin 


and charge drift-diffusion equations gnmi], 

Qij = (la) 

^ (lb) 

Jc = - e^i^CrV X (Ic) 

= 0 (Id) 

where is the charge potential (defined such that 

is the applied voltage difference), is the spin accumu¬ 

lation (the difference in spin chemical potential between 
up and down spins) for spins aligned along direction 
Qij is the spin current dyadic describing the flux along 
direction i of spins aligned along direction j, Jc is the 
charge current density, Vi is the derivative along coor¬ 
dinate direction i, eijk is the antisymmetric Levi-Civita 
tensor, a is the electrical conductivity, and Xgf is the 
spin diffusion length. These equations describe coupling 
between spin currents and charge currents such that a 
charge (spin) potential becomes a source of spin (charge) 
current. Here, we are interested in the SHE, where we 
apply a charge current to generate a spin current. In our 
analysis later in the manuscript, we ignore inverse SHE 
corrections to the charge current which appear at second 
order in Osh- The steady-state approximation used to 
obtain Eqs. is valid for time scales much longer than 
the spin relaxation time in metals (typically fs). Since the 
frequencies used in this analysis have time scales of order 
1 yus, the steady-state approximation is appropriate. On 
these long time scales, other potential carriers of angular 
momentum (spin waves, etc.) have relaxed. Einally, we 
assume that hyperfine interactions are negligible. 

The boundary conditions on the spin accumulation re¬ 
quire that no spin current flows into the surrounding in¬ 
sulator media, i.e. the perpendicular component of Q 
equals zero at the structure’s interfaces with insulating 
materials. This set of boundary conditions can be writ¬ 
ten in the form 



boundary 




boundary 


( 2 ) 


We solve Eqs. and numerically using the COMSOL 
multiphysics package, for arbitrary geometries, to gener¬ 
ate the steady-state spin accumulation and spin current 
density profiles for a given applied voltage. The net spin 
accumulation generated by the SHE is illustrated in Eig- 
ure Et for a uniform charge current flowing in a rectangu¬ 
lar thin film. The steady-state spin accumulation spirals 
around the center, increasing in magnitude toward the 
surfaces. Spins accumulate on all surfaces, polarized par¬ 
allel to the surface and perpendicular to the charge cur¬ 
rent. Due to the lossy nature of spin transport in metals, 
the spin accumulation decays exponentially away from 
each boundary within the spin diffusion length, as shown 
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in the insets of Figure [^. This figure illustrates the 
spin accumulation for spin diffusion lengths comparable 
to (blue line) and much smaller than (red line) the wire 
dimensions. Spin current, which is the spatial deriva¬ 
tive of the spin accumulation, is therefore always flowing 
in the steady-state approximation and the absence of a 
sufficiently large magnetic field. 

For a very long, rectangular geometry of thickness h, 
width re, and length L carrying a current density Jc (pro¬ 
duced by a voltage Vb over a length L), the diffusion 
equations are analytically solvable: 


cosh ( 2 ^^) 




(5) = 0 
y ^ 




Qxz 

Q ZX 


^ cosh ( 2 ^) 

esnJch I cosh(^) _ \ 

bosh (2^) ) 

OsHJch \ 

Vcosh(2i^) ) 


(3a) 

(3b) 

(3c) 

(3d) 

(3e) 


where the y-axis is defined along the length (parallel 
to the applied current), x is along the width, and z is 
along the thickness. All other components of Q are zero. 
This solution is consistent with prior results that have 
measured the interfacial spin accumulation as a function 
of thickness by its action on an adjacent ferromagnetic 
layer. The flow of spins within the material is equiva¬ 
lent to a flow of angular momentum [42]. Because total 
angular momentum is conserved, the change in angular 
momentum of the electrons must be compensated by a 
change in angular momentum of the lattice. This results 
in a mechanical torque on the material. The mechanical 
torque associated with the SHE can be calculated from 
conservation of angular momentum, assuming all spins 
gained and lost in the diffusive and SHE processes are 
converted to lattice angular momentum. Eor materials 
with no equilibrium magnetization, this assumption is 
valid. The local torque density, r, is the local net change 
of angular momentum per unit volume. Eor SHE drive 
frequencies much smaller than the spin scattering rate 
(< 10 ^^ s“^ [ISlIll]), is [45| 


Tj — ViQij (4) 

Eor a structure with dimensions much larger than the ma¬ 
terial’s spin diffusion length, the net torque, T = f f dV, 


is localized at the surfaces, 

^ 2 :, right/left S H Jc^h (5a) 

top/bottom S H Jc^/W (^t)) 

where - (- 1 -) refers to the right (left) or top (bottom) sur¬ 
face. The torque averaged over the entire structure is 
zero, so that conservation of total angular momentum is 
maintained. However, the SHE creates a separation of 
angular momentum and thus introduces local mechani¬ 
cal torques. In this limit, the total torque at each surface 
is independent of the spin diffusion length, and contains 
only one free parameter: the spin Hall angle. The torque 
at the surface is aligned parallel to the surface and per¬ 
pendicular to the applied current. This results in an ef¬ 
fective force along each surface plane, compensated by 
an opposing force distribution in the interior. 

The torque density can be converted to an effec¬ 
tive force density profile (see Supplemental Material for 
derivation) which we apply to the SHE material to cal¬ 
culate the deformation. Eor our numerical calculations, 
we use 10 nm for Xgf and 0.2 for &sh (35] |44]. Other 
properties of Pt are summarized in the Supplemental ma¬ 
terial. 

Eigure shows the deformation of a Pt film, which 
is fixed at the bottom surface, induced by the SHE. The 
SHE causes the surfaces to move in the direction parallel 
to the current while the interior moves in the opposite 
direction. Any object mechanically coupled to the SHE 
film will feel the force applied at the contacting interface. 

We simulate the mechanical effect of the SHE on a mi¬ 
croelectromechanical system (MEMS) to demonstrate its 
power as a detection method. Eqs. and are solved 
numerically using COMSOL with A^/, a, &sh^ and the 
applied voltage difference Vb as input parameters, result¬ 
ing in numerical profiles of the spin accumulation and 
spin current. The values of these parameters used in the 
simulations are summarized in the Supplemental mate¬ 
rial. 

The force profile, derived using Eq. and the pro¬ 
cess described in the Supplemental Material, is applied 
to the surface of the MEMS structure, and the resulting 
deformation is calculated using linear elastic theory in 
COMSOL. Motivated by the force profile, which is ap¬ 
proximately constant along the length of the SHE mate¬ 
rial, we choose an anchored beam, shown in Eigure [^, as 
our MEMS for this demonstration. In the future, other 
MEMS may be explored for device and detection opti¬ 
mization. The inset in Eigure illustrates the beam’s 
layered structure of Si(5,000) / Si 02 ( 1 , 000 ) / Mo(600) / 
A1N(2,000) / Mo(600) / Al 2 O 3 ( 60 ) / Pt( 20 ) (thicknesses 
in nanometers). The bulk of the beam is Si / Si 02 for the 
mechanical properties and robustness, and the beam is 
capped with a Mo/AIN/Mo piezoelectric detection layer 
[46] . The piezoelectric layer allows for precise electrical 
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(b) 




Side View 
1.85 MHz 

Side View 
4.97 MHz 



Top View 
19.18 MHz 


FIG. 2. (a) Representation of the doubly-clamped beam to 
be used in experiments (Pt and AI2O3 layers not shown). 
The inset illustrates the layers used in constructing the beam, 
(b) Several of the lowest order mode shapes for the doubly- 
clamped beam according to a numerical simulation in COM- 
SOL. 


measurements of the deformation of the MEMS, with the 
generated piezoelectric voltage being the measured out¬ 
put parameter. Other detection methods, including op¬ 
tical and magnetomotive, may in principle be employed 
instead. The thin AI 2 O 3 layer electrically separates the 
Pt from the Mo to avoid parasitic currents. The beam 
used in these simulations is 200 microns long and 40 mi¬ 
crons wide. For dynamic simulations we include Rayleigh 
damping m at a value that produces a quality factor of 
100 at 20 MHz. Figureshows representative mechani¬ 
cal eigenmodes of the beam, with eigenfrequencies in the 
rage of 10-200 MHz. A variety of eigenmodes with differ¬ 
ent symmetries exist, but only a limited number of modes 
will be preferentially excitable by the SHE, as we discuss 
later. This fact allows us to determine definitively that 
the mechanical excitation is due to the SHE. 

We first analyze the response to an applied direct cur¬ 
rent (dc). The beam reaches a final equilibrium defor¬ 
mation shown in Figure [^. In the linear elastic regime, 
the maximum deflection of the beam is approximately 
0.28 fm/mA. The output voltage is in the range of nV 
for currents small enough to avoid Joule heating (< 10 
mA), and is strongly dependent on the thickness of the 
wire and the spin flip length. Unfortunately, the thermal 


(a) Ends of Beam Fixed 



(c) 19.18 MHz 

mjjpm 56.91 MHz 
139.24 MHz 



Frequency [MHz] 


FIG. 3. (a) Deflection of the beam due to steady-state (dc) 
applied current, (b) Voltage response of oscillator when ex¬ 
cited via the spin Hall effect (blue) and piezoelectrically (red) 
when driven by a 1 mV electric potential, (c) Mode shapes of 
the three eigenmodes most strongly excited by the spin Hall 
effect. 


noise is expected to be at least an order of magnitude 
larger than the DC signal due to the SHE at room tem¬ 
perature. 

We also analyze the response under alternating current 
(ac) excitation. We apply the force profile at a given 
frequency and calculate the response at the input fre¬ 
quency, resulting in an effective scattering parameter for 
the SHE-MEMS device. Here, the force profile will only 
couple to the limited number of eigenmodes with strong 
spatial overlap. Figure If shows the response due to 
SHE excitation and due to only piezoelectric excitation. 
In both cases, the driving potential has a magnitude of 
1 mV. The piezoelectric contacts strongly excite many 
modes while the spin Hall effect only strongly excites 
three modes, at 19.18 MHz, 56.91 MHz, and 139.24 MHz. 
The three frequencies that couple strongly to the SHE 
have profiles shown in Figure [^. The two strongest peaks 
correspond to the two lowest order longitudinal modes, 
which involve vibration along the length of the wire, while 
the highest frequency peak corresponds to a more com¬ 
plicated mode shape. For the mode at 19.18 MHz, we 
can expect a response of approximately 300 ±92.5 nV 
at 300 K or 300 ±11.6 nV at 4 K. These measurements 
correspond to values for Qsh of 0.2 ±0.062 at 300 K and 
0.2 ± 0.0075 at 4 K. Since this analysis was performed 
using a 1 mV input potential and the output signal in¬ 
creases linearly with input, the signal size can easily be 
increased by at least two orders of magnitude without in- 
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creasing the noise simply by increasing the driving volt¬ 
age. The sensitivity analysis detailed in the Supplemen¬ 
tal Material. 

In conclusion, we have calculated the mechanical re¬ 
sponse of a spin Hall material to an applied current. We 
showed that the mechanical response can be coupled to 
a MEMS and, the resulting signal is large enough to be 
measured with standard equipment. The resonant behav¬ 
ior of the SHE is unique, distinguishable by the different 
modes it excites. This technique can be used to measure 
the SHE at different temperatures, for a variety of mate¬ 
rials, and can be extended to study other spin properties 
such as the surface states of topological insulators and 
magnetically-doped 3D topological insulators [48], edge 
states of quantum spin Hall insulators [49] , spin waves in 
antiferromagnets, and materials displaying the quantum 
anomalous Hall effect. Our results also offer prospects 
for new actuation methods for nanomechanical devices 
based on spin transfer effects. In-depth study of the spin 
Hall effect will reveal insights about spin-orbit coupling 
and the interplay between charge and spin transport in 
metals. 
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